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ABSTRACT 

We present results of optical spectroscopic observations of the mass donor star in SS 433 with Subaru and 
Gemini, with an aim to best constrain the mass of the compact object. Subaru/FOCAS observations were 
performed on 4 nights of October 6-8 and 10, 2007, covering the orbital phase of = 0.96-0.26. We first 
calculate cross correlation function (CCF) of these spectra with that of the reference star HD 9233 in the 
wavelength range of 4740-4840 A. This region is selected to avoid "strong" absorption lines accompanied 
with contaminating emission components, which most probably originate from the surroundings of the donor 
star, such as the wind an d gas strea m. The same analysis is applied to archive data of Gemini/GMOS taken at 
(f) = 0.84-0.30 by Hillwi£&Gies| |2008). Fr om the Subaru and Gemini CCF results, the amplitude of radial 
velocity curve of the donor star is determined to be 58.3±3.8 km s~' with a systemic velocity of 59.2±2.5 km 
s"'. Together with the radial velocity curve of the compact object, we derive the mass of the donor star and 
compact object to be Mo=12.4±l. 9 M0 and Mx=4.3±0. 6 Mq, respectively. We conclude, however, that these 
values should be taken as upper limits. From the analysis of the averaged absorption line profiles of strong 
lines (mostly ions) and weak lines (mostly neutrals) observed with Subaru, we find evidence for heating effects 
from the compact object. Using a simple model, we find that the true radial velocity amplitude of the donor star 
could be as low as 40±5 km s~' in order to produce the observed absorption-line profiles. Taking into account 
the heating of the donor star may lower the derived masses to Mq = 10.4;|;j ^ M© and Mx = 2.5;!;q^ Mq. Our 
final constraint, 1.9 M0 < Mx < 4.9 M0, indicates that the compact object in SS 433 is most likely a low mass 
black hole, although the possibility of a massive neutron star cannot be firmly excluded. 



Subject headings: accretion, accretion disks — stars: individual (SS 433, VI 343 Aquilae) 
X-rays: binaries — techniques: spectroscopic 
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1. EVITRODUCTION 

The microquasar SS 433 is a target of great interest in mod- 
ern astronomy as a unique Galactic source that sho ws steady 
relativistic (v = 0.26 c) jets (for a review see, e.g., iMargonI 
fT984t lFabrikall2004h . It gives us an ideal opportunity to 
study the formation mechanism of astrophysical jets under 
supercritical mass accretion onto a compact object. Although 
SS 433 has been studied for about 30 years since its discovery, 
the identification of the compact object, the most fundamen- 
tal issue to understand this system, has remained unsolved. In 
particular, the question whether it is a neutron star or a black 
hole is still open. 

A direct way to identify the compact object in a bi- 
nary system is to determine its mass function by measur- 
ing the Doppler shifts of the stars due to the orbital mo- 
tion. In the case of SS 433, both the inclination angle and 
the orbital period have been measured with high accuracy to 
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i = 78.8° (lMargon& Andersonlll989l) and P = 13.082 days 
("Goranskii et al.|[T998 !). respectively. Hence, if the radial ve- 
locity amplitudes of the compact object and the donor star (or 
companion star) are known, their masses can be firmly deter- 
mined. 

Most of the optical light from SS 433 is emitted from the 
compact object (i.e., from the jets and the accretion disk). 
Hence,it is relatively easy to measure the radial velocity of the 
compact object, Kx, utilizing emission lines that can be ob- 
served even with moderate-size telescopes (hereafter, the sub- 
scripts X and O represent the compact object and the donor 
star, respectively). Using He II A4686, Fabrika & Bychkoval 
(1990) obtained Kx = 175 ±20 km s"'. Their analysis is 
based on out-of-eclipse data with minimum incli nation, i.e., 
when th e disk is oriented maximally toward us. iGies et aO 
(l2002bl) used C II A723 1,7236 blended lines and determined 
Kx= 162 ±29 km s"'. The systemic velocity of the compact 
obje ct, 7x, is d ifferent between Fabrika & Bychkova (199C| 
and iGies et all (l2002b), however. IGies et al.l (l2002bl) argue 
that this is caused by the difference of the line emitting re- 
gions in the accretion disk. 

The measurement of the radial velocity of the donor star, 
Kq, is more complicated, since the signal from the donor star, 
detectable as absorption lines from its s urface, is h i dden in 
the high flux from the compact object. iGies et al] (l2002ah 
detected faint absorption lines in the blue part of the optical 
spectra of SS 433, showing Doppler shifts expected from the 
orbital motion of the donor star They found that the spec- 
trum re sembles t hat o f an A- type evolved star, as confirmed by 
ICherepashchuk et al.l (l2005h . Latest results on the radial ve - 
locity of the donor star are reported bv lHillwig & GiesI (|2008|) . 
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who used data obtained at the Kitt Peak Nati onal Observatory 
(previously published in lHillwig et alJ (120041) ') and the Gemini 
telescope. They derived a donor star semi-amplitude of Kq = 
58.2 ± 3.1 km s"' and a systemic velocity of 70 = 73 ±2 km 
s"'. Combining their results with Kx= 168 ± 18 km s~', the 
average value of Fabrika & Bychkova (1990) and iGies et aO 
(l2002bl) . they conclude the mass of the donor and compact 
object to be Mq = 12.3 ± 3.3 Mq and Mx = 4.3 ±0.8 Mq, 
respectively. In the discussion below (§ 6.1), we review in 
detail the history of the radial velocity determinations of the 
compact object and the donor star. 

The selection of absorption lines originating from the pho- 
tosphere of the donor star is a key issue for a reliable determi- 
nation of the radial velocity of the donor star It is also impor- 
tant to observe SS 433 when the disk is oriented maximally to- 
wards the observer and the outflowing m a terial d oes not inter- 
sect w ith th e Une of sight. ICharles et all (|2004|) . iBarnes et"aLl 
( I2OO6I) . and lClark et al.l (l2007h suggested that the spectral type 
of the donor star is an A-type supergiant. However, the wave- 
lengths of the lines they observed di d not follow the expected 
orbital velocity curve. iBarnes etal. (2006) and Clark etal] 
(|2007) pointed out that some absorption lines originate from 
the mass accretion flow onto the compact object, and not from 
the surface of the donor star. More importantly, the heat- 
ing of the donor surface by the compact object may signifi- 
cantly affect the accurate rn easurement of the radial velocity 
(ICherepashchuk et al.ll2005l) . 

The intensity of the absorption lines is maximum in the 
central phase of the accretion-disk eclipse and it rapidly de- 
creases when the compact object moves out from the eclipse 
([Hillwig et al. 2004). There are at least two types of absorp- 
tion lines in the SS 433 spectra, (1) absorption lines with 
emission components, and (2) pure absorption lines. The for- 
mer ones are usually stronger, showing typical "shell-type" 
line profiles, where the absorption line is located between two 
emission components. The pure absorption lines are weaker, 
and they are apparently not accompanied by emission fea- 
tures. Hereafter, we call these two types of absorption lines 
"strong" and "weak", respectively. 

The strong absorption lines have large oscillator strengths, 
and are usually formed in the upper regions of a star's atmo- 
sphere. In the case of SS 433, considering the heavy mass- 
loss rate of the donor of M - lO""* Mq yr"' (lFabrikall200"4b . 
and the underlying emission components, we expect that the 
strong absorption lines and associated emission lines may be 
formed in pla c es not directly related to the donor photosphere. 
iHillwig et al.l (|2004|) used a "highly rectified continuum" (see 
below) to smooth out the emission components of the strong 
absorption lines. Generally speaking, using the strong absorp- 
tion lines without a detailed modeling of each spectral feature, 
has to be considered as risky. 

The heating of the donor star by the strong UV radiation 
from the supercritical accretion disk is known to be important. 
For a disk UV luminosity of ^ 10"^" erg s"' , the heated surface 
of the donor star in SS 433 is expected to have a temperature 
of - 20000 K (Fabrika 2004), instead of only 9500 K in the 
absence of heating. (Hillwig e t al.l2004tlCherepashchuk et aP 
l2005h . Such a strong heating effect can produce the emis- 
sion components observed in the strong absorption lines. It 
can also distort the radial velocities measuring from the weak 
absorption lines (Cherepashchuk et al. 2005), since they are 
observed in the non-heated (or partly heated) regions of the 
donor surface, whose configuration changes with the orbital 



phase. 

Here, we present the most recent determination of the ra- 
dial velocities and mass of the donor star and the compact 
object. In our analysis, we consider various effects as dis- 
cussed above. Firstly, we carefully select absorption lines that 
originate from the donor's photosphere with minimum con- 
tamination by emission components, such as those from the 
wind, the gas stream, and the heated surface of the donor star 
For this purpose, we use high quality optical spectroscopic 
data obtained with Subaru FOCAS in 2007 October during 
four night s. Archival data taken at the Gemini telescope pub- 
lished by iHillwig & GiesI (l2008l) are analyzed as well. We 
present our best constraints on the mass of the compact object 
in SS 433, based on a model that accounts for the averaged ab- 
sorption line profiles with consideration of the heating effects 
from the compact object. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. Subaru Data 

We observed SS 4 33 with the FOCAS instrument 
dKashikawa et al.|[2002h on the Subaru telescope on October 
6-8 and 10, 2007. The jet of this source is known to pre- 
cess with a period of 162.15 days. This epoch was chosen 
to observe the system in a particular phase. Firstly, the disk 
was oriented maximally towards us « 0, where 4> is pre- 
cessional phase), which prevented the gas outflow from the 
accretion disk to intersect with the line of sight. Secondly, 
the orbital phase included the eclipse of the compact object 
by the donor star (cf) w 0, where (j) is orbital phase). The 
spectra cover the orbital phase of 0.96 < (f> < 0.26 and pre- 
cession phase of 0.02 < ip < 0.04, based on the orbital light 
cu rves presented belo w and the precession ephemeris given 
bv lGieset al.l(l2002bl) . 

FOCAS was operated with the 0".4 slit, VPH450 Grism, 
and 3x1 binning for the CCD chip. The sky condition was 
mostly photometric, except for October 10, 2007, with a typi- 
cal seeing of « 1".0. The resulting spectra cover the wave- 
length range of 3750-5250 A with a dispersion of 0.37 A 
pixel"'. Per night, we took 5 to 8 frames with 11-20 min- 
utes exposure each. As reference stars we observed HD 9233 
(spectral type A4 lab), whose spect rum is similar to the donor 
star in SS 433 (IHillwig et al.ll2004h . HD 187982 (Al la) and 
HD 332044 (B3 la). 

All spectra are reduced by the IRAF package (" Todyl [19931) 
in the usual way. We first subtract the bias, using averaged 
data of 20 bias frames. Then, to correct for the individ- 
ual difference of each frame, we further subtract the remain- 
ing offsets in the over-scanned region from the exposed re- 
gion. For the data of SS 433 with long exposures, the cos- 
mic ray particle traces are removed using the lacos_spec task 
dvan Do kkum 2001). The averaged flat image is created from 
13 frames taken each night, which is then corrected approx- 
imately for the wavelength dependence of the flux to give a 
"normalized" flat frame. Finally, we divide the object frames 
by the normalized flat frame. Any remaining cosmic ray 
traces are removed manually in this stage. 

Accurate wavelength calibration is a critical point for our 
scientific goals. We utilize a Thorium-Argon lamp with the 
identify task on IRAF. The accuracy is confirmed to be better 
than 4 km s~' by checking the interstellar absorption feature 
of Ca II A3933.66. The flux calibration is made using the 
standard stars, BDH-28D4211 (first night) and BDH-40D4032 
(second-fourth night). We ignore the effects of the slit loss, as 
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we are mainly interested in the change of the relative flux. Fi- 
nally, the atmospheric extinction is corrected. For the spectra 
of the fourth night, when the sky condition was not photomet- 
ric, we correct the flux level relative to that of the third night 
by using the B-band magnitudes reported in § 2.3. To achieve 
the best signal-to-noise ratio, we add all the individual spec- 
tra produced in this way, except for those with low statistics 
or poor observing conditions, to obtain one spectrum for each 
night. We utilize 5, 7, 5 and 7 frames for the 1st, 2nd, 3rd and 
4th night, respectively. We then produce "normalized" spec- 
tra, by dividing the original spectra by a smooth continuum. 

2.2. Gemini Data 

We analyze the archival data of SS 433 observed with 
the GMOS instrument on the Gemini telesc ope on June 7- 
13,2006 (UT). This data was also used by iHillwig & GiesI 
(j2008). The observations cover the orbital phase of </> = 0.84- 
0.30 at the precession phase of i/; = 0.02-0.06, when the ac- 
cretion disk is oriented maximally toward us. In the epoch of 
the Gemini observations, SS 433 was found to be more active 
than during the Subaru observations. Starting from the orig- 
inal frames available from the Gemini w eb site^, we red uced 
the data using the ESO-MIDAS package dWarmelsll 19921) . ac- 
cording to standard procedures. Like for the Subaru data, 
we produce an averaged spectrum for each night, "normal- 
ize" by a continuum fit, and apply a heliocentric correction to 
the wavelengths. Table 1 summarizes the observation dates 
and exposures of the Subaru and Gemini data analyzed in this 
paper. 

2.3. Photometric data 

We obtained photometric data of SS 433 in the standard B 
and V bands at the Im telescope of the Special Astrophysi- 
cal Observatory (SAO RAS) on October 2-11, 2007 with the 
CCD detector EEV 42-40. In addition, we obtained Subaru 
V-band images taken just before the spectral observation. To 
determined the B magnitudes of SS 433 during the Subaru 
spectroscopic observations, we used the V data and an inter- 
polation of the (B-V) versus V relation, wh ich is established 
very wefl in SS 433 jGoranskii et alj|1998l) . The final pho- 
tometric accuracy is 0.01 and 0.02 magnitudes in the V and 
B bands, respectively, in direct observations. For the B-band 
magnitudes interpolated to the Subaru observation time, we 
obtained an accuracy of 0.03 mag. 

Figure [T] shows the V and B photometric light curves of 
SS 433. The SS 433 brightness out of eclipse is V = 14.0, in - 
dicating that the object was in "passive state" (lFabrikall2004l) . 
The middle eclipse took place between the first and the sec- 
ond night of the Subaru observations. The minimum is very 
well shaped and regular. Using these data and all previous 
photometric data of SS 433 we have, we update the orbital 
ephemeris. The main minimum is Mini = JD 2450023.76 ± 
0.2 + (13.08227 ± 0.00008) x E. The new orbital period is 
slightly greater than the previous one (13.08211) published 
by Goranskii et al. ( 1998). This does not mean, however, that 
we detect a change of the period. The orbital ephemeris sat- 
isfy the previous photometric data as a solution with a con- 
stant period. The particular photometric eclipse displayed in 
Figure[T]has occurred 0.375 day after the predicted time from 
[Goranskii et al. (1998) and 0.18 day after the time predicted 
by our new ephemeris. Such deviations are well-known from 

' http://www4.cadc-ccda.hia-iha.nrc-cnrc.gc.caygsa/ 



the photometric behavior of SS 433. In the following, we as- 
sume JD 2454380.335 as the peak of the eclipse to calculate 
the orbital phase. For the analysis of the Gemini spectra, taken 
in June 2006, we apply the new orbital ephemeris as presented 
above. 

3. SPECTRAL FEATURES FROM THE DONOR STAR 

Figure |2] shows the flux-calibrated spectra of SS 433 in the 
3750-5250 A range obtained with the Subai'u FOCAS instru- 
ment. Apparently, the continuum fluxes were small in the 
first and second night, corresponding to orbital phases close to 
the eclipse (</) = 0.956 and 0.034), and increased as the com- 
pact object move out of the eclipse. The most prominent fea- 
tures in these spectra are emission l ines origi nating from the 
accretion disk and the gas stream (Fabrika 2004), including 
H I lines (from H/3 to Hll), He I (the strongest are 5048 A, 
5015 A, 4922 A, 4713 A, 4471 A), He H (4686 A), Fe II 
(the sti-ongest is 5169 A), and the C III/N III Bowen blend 
(« 4640 A). The broad lines produced by relativistic jet were 
not strong during our Subaru observations. The H/3" line is 
detected close to H7° line, from A 4400 A to A 4270 A 
due to the jet nutation motion. 

A large fraction of the optical emission of SS 433 origi- 
nates from the co mpact object, i .e., from the accretion disk 
and the jet bases ( lFabrikall2004l) . When the accretion disk 
move out from the eclipse, absorption lines from the donor 
star become very weak and hence careful analysis is required 
to study their features. To measure the orbital motion of the 
donor star (§ |4]l, we need to determine the cross-correlation 
function (CCF) with the spectrum of a reference star first. 
For this, we define three different spectral regions that are 
not affected by prominent emission lines from compact ob- 
ject; Region 1 (4490^630 A), Region 2 (4740^840 A), and 
Region 3 (4950-4990 A). 

The normalized spectra of Region 1 and Regions 2-3, taken 
during first night, are plotted in Figures [3] and |4] respectively. 
Region 1 contains many strong absorption lines of Fe II and 
Ti II surrounded by emission components. By contrast. Re- 
gion 2 is practically void of strong absorption lines and con- 
tains the Cr II A4824 line with a weak emission component. 
Region 3, the narrowest one, contains the Fe I A4957 line with 
a weak emission component. To make the absorption fea- 
tures clearly visible, we further divide the "normalized spec- 
tra" by a continuum function modeled by Legendre polynomi- 
als of order «15 in each region. We call the resultant spectra 
"highly rectified spectra". 

Figures |5] and |6] show the highly rectified spectra of SS 433 
in Region 1 and Region 2-3, respectively, together with the 
normalized spectrum of the standard star HD 9233. HD 9233, 
spectral type A4 lab, is known to show an a bsorption line 
spectrum similar to the donor star of SS 433 (iHillwig et al.l 
2004). For an easy comparison, all the spectra have been 
shifted into the rest frame by correcting for its radial veloc- 
ity as determined by the CCF analysis in the next section. It 
is known that even during the eclipse the surroundings of the 
compact object (probably the accretion disk wind) contribute 
significantly to the total brightness of the system. The spec- 
ti-um of HD 9233 is scaled to match the flux of flie SS 433 
spectra by multiplying with a factor of 0.36 (IHillwig et al.l 
120041) . The deep absorption lines at 4500 A, 4 760 A, 4780 
Aand 4980 A are due to interstellar absorptions (iHobbs et al.l 
120081) . Apart from these, the specti-a of HD 9233 and SS 433 
contain the same set of absorption lines. These absorption 
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features in the SS 433 spectra bec ome deeper a s the donor 
star hides the compact object (Hill wig et al.ll2004l) . providing 
evidence that they originate from the donor star. 

The Gemini spectra observed for 7 nights are analyzed with 
the same procedure. Figure [T] shows the resuhing highly rec- 
tified spectra in Regions 2 and 3 obtained with the Gemini 
GMOS. For comparison, we also plot the Subaru HD 9233 
spectrum in the same figure. The wavelengths are corrected 
for its radial velocity, except for that of the seventh night 
where the absorption-line features are found to be extremely 
faint. We do not analyze the Gemini spectra of Region 1 , since 
the strong absorption lines in the region appear as pure emis- 
sion lines. This is probably due to the higher activity during 
the Gemini observations than during the Subaru observations. 
This make it impossible to use them for our CCF analysis. 

4. THE RADIAL VELOCITY OF THE DONOR STAR 
4.1. Cross Correlation Function Analysis 

We derive the radial velocity of the donor star by cross- 
correlating the spectra of SS 433 with those of HD 9233 in 
each spectral region (Region 1, 2 or 3). Thereby, we assume 
the heating effects by the compact object at the surface of the 
donor star are negligible (see discussion in § 14.31 ). We further 
ignore wavelengths with strong interstellar absorption. We 
pay special attention to the determination of the radial velocity 
of H D 9233, which is known to be a radial velocity variable 
star (iHiUwig & GiesI 2008). By measuri ng the Doppler shifts 
of 13 non-blended absorption lines (Chentsov & SarkisvanI 
l2007h fr om the A4 lab supergiant, we derive a systemic ve- 
locity of HD 9233 of 7hd 9233 = -44.2 ± 1 .3 km s"' . To verify 
our analysis, we also measure the radial velocity of another 
reference star HD 187982 (Type Al la), which was observed 
on October 7, 2007 (i.e., during the second night of our Sub- 
aru observation s). Our result is in good agreement with the 
literature value (IWilsonlll953h within the error range. 

Figure |8] shows the radial velocity curve of the donor star 
in SS 433 from the CCF analysis of the Subaru spectra. The 
error bars in each point correspond to the "standard error" 
of the CCF (lFernielll989l) . The amplitudes of the radial ve- 
locity curve are different between the three spectral regions 
used in our analysis, with Region 2 showing the largest am- 
plitude and Region 1 showing the smallest one. In Figure |9l 
we display the CCF results for the Gemini spectra of SS 433 
together with the Subaru spectrum of HD 9233. We confirm 
that the amplitude of the radial velocity curve derived from 
the Gemini spectra shows the same tendency as in the Subaru 
case (i.e., the amplitude derived from Region 2 is greater than 
the on from Region 3). Table |2] summarized our results on the 
radial velocities. 

We have demonstrated that the selection of absorption lines 
does affect the estimate of the radial velocity within our sim- 
ple analysis. The differences in the amplitudes are related 
to the strength of the spectral features; the absorption lines 
in Region 1 are the deepest and have underlying emission 
components, while those in Region 2 are the weakest and 
mainly do not show emission components. We interpret that 
the strong absorption lines are more significantly affected by 
the emission from the wind, the gas stream, and the heated 
surface of the donor star, which decrease the amplitude of the 
radial velocity curve. From Region 1, we obtain = 24 ± 9 
km s"' with a systemic velocity of 70 = 52 ± 6 km s~' . This 
value for Kq is even smaller than the result from lHillwig et al.l 
(I2OO4) . = 45 ±6 km s"', who studied the same spectral 
region. This is probably due to different conditions of the sur- 



roundings between the two epochs of observations. 

In this context, the selection of "weak" lines is important to 
determine correctly the motion of the donor star, as for as the 
region responsible for the production of the absorption lines 
is constant over the orbital phase. Under this assumption, we 
can estimate the amplitude of the radial velocity curve by fit- 
ting the velocities with a Keplerian solution. FigurefTOlshows 
the Subaru and Gemini results obtained from the CCF analy- 
sis of Region 2 together with the best-fit curve. We obtain a 
semi-amplitude of the radial velocity of Kq = 58.3 ± 3.8 km 
s"' and a systemic velocity of 70 = 59.2 ±2.5 km s~'. This 
value of Kq is consistent with the result of Hillwig & GiesI 
(120081) within the error bars. We note that Figure [TO]may in- 
dicate a distortion of the donor's radial velocity curve in the 
orbital phases 0.0-0.15. 

4.2. Average Absorption Lines Profiles 

The high quality of the Subaru spectra allows us to stu dy in - 
dividual absorption lines. In the case of CCF analysis (§ 14. lb . 
one compares two stars (i.e., of SS 433 and reference stars) 
with almost identical spectra. Hence, the complex blending 
and crowding of the absorption lines is not critical for the 
study. The analysis of the individual lines, however, depends 
strongly on such line blending effects and requires knowledge 
of the "laboratory" wavelengths of the blends. We carefully 
check the whole spectrum of SS 433 for the first and second 
night, when the system was in maximum eclipse. First, we 
select two groups of absorption lines: "strong" lines, which 
clearly show emission components, and "weak" lines with 
pure absorption line profiles. For the line identification, labo- 
ratory wavelengths and relative line strengths, we refer to the 
Atomic Spectra Database^ and to the Atomic Line List**. We 
estimate the effective wavelengths of the blends by weighing 
the wavelengths of the individual lines according to their line 
strengths. We test each line or line blend to have the same ra- 
dial velocity in the given Subaru night, allowing a difference 
of up to 20 km s"' . The set of strong lines as well as the set of 
weak lines are both formed from lines that have same radial 
velocities within the first and the second night. The data of 
the other nights are not considered at the line selection. Ob- 
vious or resolved line blends are not included in the two line 
groups. 

Finally, we add the line profiles for each group in the 
normalized spectra in the radial-velocity space to create the 
strong and weak average line profiles. In the averaging pro- 
cedure, we apply a weight of unity for a single line and a 
smaller weight for obvious (but non-resolved) blends. The 
line blending may distort the line profiles. The averaging pro- 
cedure minimizes this distortion, because the line blending is 
only occasional, and because obvious and strong blends are 
excluded from the set of lines. 

To derive the average line profiles, we do not use the 
"highly rectified" spectra. Instead, we apply a linear contin- 
uum rectification to the final average line profiles in order to 
subtract the continuum levels near the lines. Both the "strong" 
and the "weak" absorption lines are considerably fainter than 
the strong emission lines in the SS 433 spectrum. Hence, in 
case where an absorption line is located near the wing of a 
strong and broad emission line, its local continuum becomes 

' National Institute of Standards and Technology; 
http://physics.nist.gov/PhysRetData/ASD/ 

" Department of Physics and A stronomy, University of Kentucky; 
|http://www.pa.uky.edu/~ peter/atomic/ 
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not flat. Applying a linear function for the continuum rectifi- 
cation of the average spectra makes it possible to obtain a flat 
normalized continuum for the average line profiles. Using our 
B-band photometry results (Figure [T]i, we scale the average 
line profiles to a relative intensity unit, where the contribu- 
tion from the donor star can be directly compared for the four 
nights. In the scaling process, we multiply the average pro- 
files by coefficients depending on the system brightness and 
adopt the first Subaru night coefficient as 1 .0. 

We include 8 individual lines in the average strong line 
profile (Mg II A4481.21, Ti IlH-Fe II A4549.63, Ti IH-Fe II 
A4555.49, Fe II A4576.39, Ti IH-Fe II A4583.41, P IH-Cr II 
A4823.84, Si II A5041.03, Fe I A5226.86) and 8 individual 
lines in the average weak line profile (Cr I A4 16 1.42, Fe I 
A4271.76, Ti II A4290.22, Ti I A4325.13, Fe I A4528.87, 
A4983.85, A5125.il, Mg I A5183.60). We find that nearly all 
strong absorption lines with emission components are gener- 
ated by ions, while most of the weak, pure absorption lines are 
produced by neutral atoms. This implies that the ion absorp- 
tion lines are partial formed in a more extended gas envelope 
of the donor star (i.e., the donor's wind). This wind could be 
a low velocity wind, since the donor overfills its Roche lobe. 
The emission components in the strong absorption lines may 
be partly formed in the gas stream (probably the mass flow, 
see below), which is best observed in hydrogen and He I emis- 
sion lines (Crampton & Hutchings 1981; Fabrika et al. 1997; 
i:abrika 2004). 

The average strong and weak line profiles are shown in the 
right hand panels of Figures [TT| and [T2l respectively. The 
strong lines show clear evolution of their emission compo- 
nents over the four nights of observations. The absorption 
components shift with time to positive velocities with an am- 
plitude of « 40 km s"' . The weak lines show the same system- 
atic evolution with orbital phase, although its orbital shift in 
radial velocity is notably larger than that of the strong lines. 
The final signal-to-noise ratio in the average line profiles is 
very high. We detect several features in the line profiles that 
change from night to night in ways that are difficult to inter- 
pret. They would be due to line-emitting and -absorbing re- 
gions in the system having a complex structure. Additionally, 
emission components could partially fill the absorption pro- 
files of the weak lines as well. In the following, we study the 
main features only, such as the line positions and intensities. 

We measure the absorption line positions for the average 
weak and strong line profiles and compare them with the CCF 
results. Note that the CCF analysis and the study of the aver- 
age absorption line profiles are independent. We find that the 
weak absorption lines show the same behavior as the lines in 
the CCF analysis of Region 2, which is free from strong lines. 
The total radial velocity amplitude measured between the first 
and 4th spectrum is w 63 km s"', which is very close to w 67 
km s"' derived from the CCF analysis (Table |2]i. The abso- 
lute velocities are about the same as well, with the difference 
between the two methods being within 5 km s"' . We further 
compare the behavior of the strong absorption lines with the 
CCF results of Region 1 , which mainly contains strong lines 
with emission components. The difference is 10 km s"' for 
the total radial velocity amplitude, although the systemic ve- 
locity is larger by 15 km s"' for the average strong line profile 
than for the CCF analysis. Note that for the CCF analysis, 
we used the "highly rectified" spectra, where smoothing of 
the emission components may produce a systematic shift in 
the absorption line position. In the following, we study the 
average line profiles using a simple model of a close binary 



by taking into account the heating effects from the compact 
object. 

4.3. Study with a Simple Model Including Heating Effects by 
the Compact Object 

4.3.1. Model Description 

Figure [13] sketches the binary system with its main compo- 
nents. Since the UV lumi nosity of the compact object is as 
high as Luv ^ 10"*" erg s~' (lFabrikal2 004). the compact object 
can heat the donor surface up to ~ 20000 K during the 13- 
day binary period. Although the detailed of geometry of the 
system is unknown, we know that the size of the optical con- 
tinuum source is about the size of the donor-star or slightly 
exceed it, because optical eclipses are never total (about half 
of the continuum light remains present always). The donor 
star has an extended and dense envelope, which is probably 
due to a strong, low-velocity wind. Studies of X-ray eclipses 
in SS 433 (Filippova et al. 2006) revealed that the radius of the 
envelope that is opaque to X-rays exceeds the donor radius (or 
its Roche lobe radius) by 10-20 %. This proves the existence 
of a gas envelope ("coat") around the donor star, which can 
produce the emission components observed in the strong ab- 
sorption line profiles. This gas envelope is also sketched in 
the figure. 

We construct a simple toy model of the system to study 
the heating effect on the absorption lines. Heating effects in- 
crease flie_observedradial_veloci^^ absorption 
lines jAntokhina et alj|2005l; ICherepashchuk et al.1 2005). be- 
cause, due to the spin of the donor star, the non-heated side 
moves with a larger velocity than that of the center of mass. 
In out model, we consider three different regions of the donor 
surface (Figure [T3]i: Region I, which is not heated and pro- 
duces absorption lines only. Region II, which is heated and 
produces emission lines only, and Region III, which is over- 
heated and therefore does not produce any emission or absorp- 
tion lines of elements with low ionization potentials (like, i.e., 
Ti II or Fe II). 

We assume that the system is synchronized and the donor 
star is a sphere with a volume equal to that of the Roche lobe 
for a given mass ratio q = Mx/ Mq. The orbital inclination 
of the system is adopted as 79° ( iFabrikall2004l) . The semi- 
amplitude of the compact object is set to Kx = 160 km s"', 
as derived late in this paper ($|5]l. The semi-amplitude of the 
donor star Kq (or the mass ratio q) is taken as free parameter 
in this model. The donor surface is divided into 100 grid cells 
both in longitude and latitude, with each grid cell producing 
a Gaussian absorption line profile in the non-heated region 
or an emission line profile in the heated region. A Gaussian 
line width of FWHM = 5 km s"' is adopted. For absorption 
lines formed at the donor's surface which is not exposed to 
the UV source (Region I in Figure [T3]l. the Gaussian line in- 
tensity is normalized to unity. For emission lines formed in 
the heated region of the surface (Region II in Figure [T3]l. we 
calculate the normalization depending on the heating param- 
eters (see below). We ad opt a quadratic limb-darkening law 
jKallrath & Milon3ll999h with x = y =\ for absorption lines 
only. Using this limb-darkening law, we can easily fit the 
observed absorption lines. Note, however, that determining 
the correct limb-darkening law for this supergiant with heavy 
mass loss and heating is a truly complex task. In the case of 
emission lines, we do not account for a limb-darkening ef- 
fects, since one expects an inverse temperature gradient in the 
emission region. 
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The gas envelope is modeled in a similar way than the donor 
surface, with the difference that it produces emission lines in 
the heated region (Region II), but no emission or absorption 
lines in the other two regions. The radial extent is set to 10 % 
of the donor radius. It has 10 individual segments in the radial 
direction. We assume that the emitting gas in the envelope 
rotates with the same velocity as the donor surface and moves 
in radial direction with the escape velocity Vesc. 

The UV source is spherical and has the same size as the 
donor star in the model (Figure [T3]). Each point of the donor 
surface sees all the points of the extended source visible from 
it. We calculate the angle of incidence of the UV radiation 
in each point of the donor surface. Thus, we specify only 
geometrical properties of the source. We suggest that the 
donor's regions exposed to the UV radiation of the source pro- 
duce emission lines. This is expected if the temperature gra- 
dient i n the donor's atmosphere is inverse. lAntokhina et alj 
( l2005h confirmed this behavior in their X-ray heating model 
of a low-mass X-ray binary. UV radiation is subject to strong 
extinction and may therefore not penetrate deeply into the 
donor atmosphere. Considering the high luminosity (Luv ^ 
10'*°erg/s) of the accretion disk in SS 433, however, we sug- 
gest that the UV radiation can indeed produce the inverse tem- 
perature gradient in the atmosphere and that the gas in the 
donor's wind may be ionized down to the photosphere. 

We introduce two empirical parameters, the heating coeffi- 
cients 7phot and 7env, which define the intensities of the emis- 
sion lines formed in the heated regions of the donor surface 
and in the envelope, respectively. From the relative fluxes of 
the incident radiation in each point exposed to the UV source, 
we calculate the expected emission-line components. We then 
determine the values of 7phot and 7env by comparing the line 
profiles between the model and the data observed in all four 
nights. The values are independent for the strong and the 
weak lines (in particular, 7env = for the weak lines). The 
obtained heating coefficients are relative ones only, and can- 
not be used to estimate the heating effects physically. To pro- 
duce a reasonable agreement between the observed and the 
modeled line profiles, 7phot has to be 2-3 times larger for the 
weak lines than for the strong lines. The weak lines do not 
require an additional emission component formed in the en- 
velope. In case of the strong lines, this component (7env) is 
necessary, since it produces emission line wings which are 
notably broader than the photospheric line profiles. By vary- 
ing other parameters of the model, such as the mass ratio, we 
can infer the required amount of heating that is necessary to 
account for the observed line profiles. 

Finally, we integrate the line profiles from individual re- 
gions of the donor surface and of the envelope which are vis- 
ible from the observer during the orbital phases of the Subaru 
observations (indicated in Figure fljTl. The final line profiles 
are convolved with the instrumental response, which is de- 
rived from single line measurements in the comparison-lamp 
spectra. The normalization of the final absorption lines is de- 
termined from the data of the first night only. These normal- 
ization coefficients are kept constant for all four nights of the 
observations. 

The final line profiles are notably broader (FWHM 
100 km s"') than our Gaussian lines formed in the individ- 
ual surface cells and they do not depend on the adopted line 
width of the individual lines (when it is less than 20km/s). 
This simple approach is sufficient for the following study, be- 
cause we do not compare in detail the observed and calculated 
line profiles, but investigate the main features of the heating 



effects only. 

Several effects are not taken into account here. In reality, 
heating by the UV radiation is a complex process. For exam- 
ple, the UV photons may not reach the donor surface because 
of strong absorption in the wind. The UV absorption can heat 
the gas deeply down to the surface, however. We further as- 
sume isotropy of the UV source, which is probably are over- 
simp lification. F or instance, the thick outer rim of the disk 
(Filip pova et al.l (2006) may cast an extended shadow. This 
effect is most important in the fourth Subaru night, since, at 
precession phases ip ^Q, the donor star crosses the disk plane 
at (j) «0.25 and 0.75. The assumption of isotropic UV ra- 
diation will lead to an overestimate of the heating effects in 
this case. Since the environment of the compact object (i.e., 
the wind, the jet bases and the disk structure) are basically 
unknown, a complex modeling of the heating is problematic. 
Thus, in this paper, we illustrate how the heating distorts the 
radial velocities of the absorption lines in SS 433 for a better 
understanding of the principal difference between the strong 
and weak lines. 

The gas strea m is a strong source of hydro gen and He I 
emission lines. ICrampton & Hutchin^sl (Il98ll) showed that 
the radial velocity curves have the largest redshifts at or- 
bital phases ~ close to the inferior conjunction of the 
donor star; their orbital phases lag behind the accretion- 
disk phase by 0.2- 0.25. Later stu dies of the hydrogen and 
He I emission lines jKopylov et al.l il989: Fab rika et al .11 19971: 
iGoranskii et al.lll997h revealed that H/3 and He I radial veloc- 
ities show the largest redshifts at orbital phases 0.9-0.95 and 
0.85-0.9, respectively. These allows also detected a partial 
eclipse in the hydrogen emission lines at orbital phases 0.1- 
0.2. This indicates that the He I and hydrogen emission lines 
in the SS 433 spectra are formed in the gas stream onto the 
a ccretion disk. 

IGoranskii etal] ( Il997h and iGies et alj (l2002bl) also dis- 
cussed whether the behavior of the hydrogen and He I ra- 
dial velocity curves may result from an evacuation of the disk 
wind surrounding the donor star, which leads to anisotropic 
wind and the observed radial velocity curves. This cannot ex- 
plain, however, both the partial eclipses and the differences 
between the hydrogen and He I radial velocity curves. If hy- 
drogen and He I lines are formed in the anisotropic wind, their 
radial velocity amplitudes have to be greater than that of the 
accretion disk (traced by He II line), since the accretion disk 
powers the wind. This does not agree with the observed radial 
velocity amplitudes. We therefore conclude that the hydrogen 
and He I emission lines are formed in the gas stream, although 
a fraction of this emission may be formed in the disk wind as 
well. In any case, the emission region must be extended. A 
probable location of the gas stream region is shown in Fig- 
ure [13] It is noteworthy that a fraction of the emission of 
the strong absorption lines (Figures [TT]) may be formed at the 
same location as the hydrogen and He I lines. 

During all for nights, we detected an additional red emis- 
sion line component in the strong and even in the weak ab- 
sorption lines, which we fail to reproduced with our model. 
We ascribe this component to the gas stream, which con- 
tributes stronger to the red emission in the strong lines than 
in the week lines. We do not model any probable eclipses 
of the gas s tream (" Kopvlov et all 1 19891: iFabrika et all 1 19971: 
iGoranskiie t al. 1997), which may change the intensity of the 
red emission components. To model this additional red emis- 
sion, we decide to follow the ra dial velocity curves o f hydro- 
gen and He I emission lines by iFabrika et al.l (Il997h . where 



Optical Spectroscopy of SS433 



7 



the lines show the largest redshifts (100-150 km s~') at the or- 
bital phases ^0.85-0.95 and the velocity decreases gradually 
with the orbital phase. Finally, the gas stream is modeled to 
produce Gaussian emission lines, whose parameters are tuned 
to reproduce the observed line profiles. 

4.3.2. Comparison with the Data 

The left hand panels of Figures [TT] and [12] show the best-fit 
models of the strong and weak line profiles, respectively, to be 
compared with the observed ones in the right panels. These 
line profiles are modeled using a radial velocity amplitude of 
the compact object of Kx= 160 km s"' and the real radial 
velocity amplitude of the donor star Kq = 40 km s"' (i.e., q = 
0.25). The relative heating coefficient 7phot for the weak lines 
is twice as large than that of the strong lines, and the wind 
velocity is set to Vesc = 260 km s"' (for the strong lines only). 
We see that our model can reproduce the overall features of 
the observations. For the strong absorption lines, the emission 
components evolve in agreement with the idea that the gas 
envelope, which rotates with the same velocity as the donor 
star, is heated by the compact object. 

For the weak absorption lines, emission lines originating 
from the heated donor surface are required, while those from 
the envelope are not. These emission lines are seen as a low 
intensity emission bump near the absorption line (as for the 
strong lines). This emission component alters the position and 
intensity of the absorption lines. 

The absorption lines move in accordance with the orbital 
phase, and their velocity amplitudes are smaller for the strong 
lines than for the weak lines. The absorption line intensities 
generally decrease with orbital phase because of the heating 
(note that the observed spectra are scaled from the photomet- 
ric data in order to keep the non-illuminated continuum ra- 
diation of the donor constant). Naturally, the observed ab- 
sorption line profiles are more complex than those produced 
by our simple model. For example, the weak absorption line 
profile of the first night (solid line in Figure [T2]| shows either 
an additional absorption in its red wing or an additional emis- 
sion component. Such a feature is not produced in our model, 
since we do not take into account a possible absorption of the 
continuum radiation from the compact object in the donor's 
wind. This effect might be important during the first and sec- 
ond night (see Figure[T3]l. Additional absorption features may 
be present for the second and third night in the blue wing of 
the strong and weak absorption lines (Figures fTTl fT2l i. Dur- 
ing these orbital phases, the speculated wind from the donor, 
which is seen projected onto the strong continuum source, is 
directed toward us because of the stellar rotation. This fea- 
ture may also produce the distortion of radial velocity curve 
observed in Figure [TOl 

We conclude that the model reproduces both the intensity 
and the radial velocity variations of the emission components 
and the absorption lines. The modeled line profiles are in- 
fluenced mostly by the following parameters; the real veloc- 
ity amplitude of the donor star Kq, and the heating efficiency 
coefficients 7phot and 7env As mentioned above, the best-fit 
value is Kq k, 40 km s"' (i.e., q « 0.25 with Kx ~ 160 km 
s"'). This is 18 km s"' or 30 % less than that measured in our 
CCF analysis (58.3±3.8 km s"'). It is impossible to produce 
the average line profiles in our model for Kq < 35 km s"' . 
Further, for velocity amplitude of the compact object between 
170 and 150 km s~', we obtain values of Kq between 35 and 
45 km s"'. We thus conclude that the donor's real radial ve- 
locity amplitude is = 40±5 km s"', based on our simple 



model. 

From the set of models reproducing the observed averaged 
line profiles, we find that the size of the overheated region 
(Region III in Figure [T3] l, modeled as a cone with origin at 
angle of the donor, has a half-opening angle of « 15-20° for 
the strong lines and «20-25° for the weak lines, respectively. 
The wind velocity of the donor is Vesc 260 km s"' . The mass 
flow, which produces parts of the red emission components in 
the absorption lines, has a radial velocity decreasing from 160 
km s~' in the first night to 90 km s~' in the last night for the 
strong absorption lines, and from 100 km s"' to 80 km s"' 
for the weak also lines, respectively. The emission line com- 
ponents formed in the gas stream are broad with a FWHM 
of 140-200 km s"'. The intensity of this component is two 
times weaker in the weak absorption lines than in the strong 
absorption lines. The origin and formation of the gas stream 
(mainly observed in the hydrogen a nd He I emission lines) 
were discussed in previous pa pers fcr ampton & Hutchingsl 
I198U iFabrikaet al.lll997l: l^rika 2004). It is important to 
note that the introduction of the red emission components in 
the modeled line profiles is necessary, although its formation 
remains unclear All these parameters do not change, how- 
ever, the modeled absorption line profiles so strongly as the 
real radial velocity amplitude of the donor and heating effi- 
ciency coefficients do. 

Finally, in Figure [141 we present the radial velocity curves 
derived from the observed average line profiles, simply based 
on the position of the absorption line minimum. For compar- 
ison, the values from our best-fit models (Kq = 40 km s~' and 
Kx = 160 km s"') are displayed as well. Clearly, the heating 
model reproduces well the observed radial velocities. Note 
that these velocities are apparent ones and differ from the real 
radial velocities of the donor star considered in the model. 

5. RADIAL VELOCITY OF THE COMPACT OBJECT 

To constrain the radial velocity of the compact object, we 
analyze the He II A4686 emission line, the brightest line 
known to originate from the compact object. Figure[T5]shows 
the line corresponding profiles. The flux level is normalized to 
that of the first night (0=0.956), calibrated using the B-band 
magnitudes (Figure [TJ. Although the line profile is complex, 
it is obvious that the line center moves from the red to the blue 
with increasing orbital phase. 

To determine the radial velocity of the He II line, we calcu- 
late the center of gravity above a certain threshold in order to 
discard the broad wings. The line wings are stronger in the 
red than in the blue. We estimate the error of the so-derived 
radial velocities by changing the flux thresholds (upper and 
lower) used in the calculation of the center of gravity. The 
results are summarized in Table [3] 

In Figure [16] we show the radial velocities of the compact 
object as measured from the He II line. A large velocity of 
w 150 km s"' is required to fit them with a Keplerian ve- 
locity curve. This is unlikely and probably due to the fact 
that the He II line was significantly affected by the eclipse 
during the first three ni ghts, as observed in previous studies 
jFabrika & Bvchkova '1990). Indeed, the effects of the eclipse 
are clearly seen in Figure [15] where the He II line profile 
changed notably across the eclipse of the line emitting region 
by the donor star If we use the data of the fourth night only, 
together with a fixed systemic velocity of 70 = 59.2 km s~', 
the same value of 7x for the donor star, we obtain Kx=159±7 
k m s~' . This value is consisten t with previous results. 

iFabrika & Bvchkoval d 19901) reported = 1 75 ± 20 km s"' , 
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using the He II A4686 line observed in th e precession phase of 
0.9 < V' < 0. 1, but outside of the eclipse. iFabrika et all (Il997h 
constructed precessional and orbital radial velocity curves of 
the He II line using additional spectral data. They found Kx = 
176 ± 15 k m s~' for th e same precessional phase of 0.9 < V' < 
0.1, while Gies et al.i (.2002b.) used C II A723 1 , 7236 blended 
lines and derived Kx = 162 ± 29 km s"' . For consistency, we 
adopt ^rx=168±10 km s~' in this paper, the average between 
these three studies a nd o ur own estimate of the K x ■ Note that 
iHillwig et al] (|2004 and iHillwig & GiesI ( |2008|) adopted the 
sa me value, 168 ±18 k m s~ ' , as the average between two stMd- 
ies. lGies etan(l2002bh and lFabrika & Bychkoval(ll990h . 

6. DISCUSSION 

6. 1 . Review of the Dynamical Mass Determination ofSS 433 

In this subsection, we review recent work on the dynamical 
determination of the mass function of SS 433 from measure- 
ments of Kx and Kq, which we compare with our results. 

1 . iGies et akl (l2002bh interpreted that an absorption fea- 
ture in the strong emission line of He I A6678 may orig- 
inate from the donor star, using the spectra taken at the 
KPNO 0.9 m telescope. They derived Kq = 126±26 
km s~' . Combining this values with Kx=n5±2Q 
km s"^ (iFabrika & Bychkovalll990l) . they determined 
Mx=16±6Mq andMo=23±8 Mq. 

2. iGies et al] (l2002ah measured the radial velocity of 
the donor star using the CCF technique for the first 
time, based on data taken with the 2.7 m telescope 
of the University of Texas McDonald. For the 4060- 
4750 A range, they calculated CCFs between different 
SS 433 spectra. They obtained A'o=100±15 km s"' 
and 7o= -44±9 km s"'. Together with K^=n5±20 
km s~' (IFabrika & Bychkoval 119901) . they concluded 
Mx=ll±5 Mq andMo=19±7M0. 

3. IHillwig et all (|2004 appUed a CCF analysis of Re- 
gion 1 (as defined in our paper) to spectra taken at 
the KPNO 4 m telescope (in this epoch. Region 1 did 
not contain strong em ission lines, unlike in the case of 
IHillwig & Giesll2008l) . They used the same reference 
star HD 9233 (type A4 lab) as we used in this paper, 
and derived A'o=45±6 km s"' and 7o=65±3 km s~'. 
Together with /rx=168±18 km s~' (the average value of 
IFabrika & Bvchkova .l99a and lGies et a l. 26o2b), they 
determined Mx=2.9±0.7 Mq and Mo=10.9d=3.1MQ. 

4. ICherepashchuk et"an (l2005l) examined Doppler shifts 
of 22 absorption lines in the 4200-5300 A range rel- 
ative to the laboratory frame, using spectra taken at 
the SAO 6 m telescope. They obtained Ab=132±9 km 
s"' and 7o=14 k m s~'. Combining these v alues with 
K^=n5 km s-i (IFabrika & Bvchkova|[T990l) . they de- 
termined Mx ~ 18 Mq and Mq «24 Mq. They noted 
that this radial velocity semi-amplitude is probably an 
upper limit, because the heating of the donor star in- 
creases the observed amplitude. 

5. iBarnes et alj (120061) observed SS 433 with the Calar 
Alto Observatory 3.5 m telescope, and the Observatory 
del Roque de Los Muchachos 2.5 m and 4.2 m tele- 
scope. They cross-correlated the SS 433 spectra in the 
4500-4630 A range with those of HD 9233, using the 



radial velocity given in IHillwig et alJ ([200l . They ob- 
tained /:o=69±4 km s"' and 7o=-53±3 km s"'. But, 
their observations were performed when the accretion 
disk was close to an edge-on orientation, where the out- 
flowing material produces strong absorption lines. Note 
that they assumed a ve locity for HD 9233 of -34 km 
s"' (iHiUwig et al.l2004l) . Since this star is velocity vari- 
able, its velocity might have been different at the time 
of the observations, and the systemic velocity 70 has an 
additional factor of uncertainty. 

6. IHillwig & GiesI (|2008|) made a CCF analysis in Re- 
gions 2 and 3 using the Gemini data. They ob- 
tained Ko=5^.2zL3.1 km s"' and 7o=73±2 km s~', 
which tog ether with.gx=16 8ib l8km s~' (the average of 
IFabrika & Bvchkovall 19901 and lGies et al. 2002b) leads 
toMx=4.3±0.8MQ andMo=12.3±3.3MQ. 

6.2. Constraints on the Mass of the Compact Object 

Once the amplitudes of the radial velocity curves of both the 
donor star and the compact object are known, one can deduce 
the mass of each component. By fitting the radial velocity 
curve of the donor star obtained from the CCF analysis with a 
Keplerian curve (i.e., without consideration for the heating ef- 
fects), we derived Kq = 58 .3 ±3.8 km s~', which is consistent 
with the value obtained byl Hillwi g & GiesI (12008). This is not 
surprising, since we use the same Gemini data for our anal- 
ysis in addition to the Subaru data, although we restrict the 
wavelength range to Region 2 only, in order to avoid system- 
atic effects from emission components in the absorption line 
profiles. We adopt the amplitude of the radial velocity of the 
compact object of /rx=168±10 km s~' and conclude the mass 
of the donor star and compact object to be Mq = 12.4 ± 1.9 
Mq and Mx = 4.3 ± 0.6 Mq, respectively. The corresponding 
mass ratio is ^ = 0.35. Again, these values should be taken 
as upper limits only if we consider the heating effect, as dis- 
cussed in Section l43] 

By taking into account the heating effects, we derive lower 
limits on the masses, since we assumed a real radial velocity 
amplitude of the donor star of = 40 ± 5 km s"' . This leads 
to lower masses of Mq = 10. 4+^ 9 Mq and Mx = 2.5:l:j}g Mq. 
We thus conclude that the compact object in SS 433 is most 
likely a low mass black hole. However, the possibility of a 
massive neutron star cannot be firmly ruled out at present, 
given the fact that a neutron stars could have masses of up 
to 3 M q as inferred from both theory ( Latt imer & PrakashI 
12007!) and observations iFreire et all2()08l) . 

7. CONCLUSIONS 

1. To study the radial velocity curve of the mass donor 
star in SS 433, we obtained high quality optical spec- 
tra with Subaru/FOCAS, covering the orbital phase of 
(j> = 0.96-0.26. We combine the se observa t ions w ith 
the Gemini data reported by Hifl wig & GiesI (l2008l) to 
analyze the largest set of the best quality spectra ob- 
served right now from this source. This allows us to 
study in detail the behavior of the "weak" absorption 
lines from the donor surface, which are least affected 
by the emission components from the surroundings of 
the donor star. 

2. We demonstrate that the selection of the spectral re- 
gion is critical for the cross correlation function (CCF) 
analysis. We adopt the 4740^840 A range (Region 2) 



Optical Spectroscopy of SS433 



9 



for this study, where only "weak" absorption lines from 
the surface of the donor star are present. If we instead 
use the 4490^630 A range (Region 1), which contains 
many "strong" absorption Unes associated with emis- 
sion components, we obtain a significantly smaller ve- 
locity amplitude than Region 2. 

3. From the Subaru and Gemini CCF results (Region 2), 
we determine the amplitude of the radial velocity curve 
of the donor star to be 58.3±3.8 km s"' . Together with 
the radial velocity of the compact object, 168±10 km 
s~', we derive masses of the donor star and the com- 
pact object ofMo= 12.4 ±1.9 Mq and Mx = 4.3 ±0.6 
Mq, respectively. These values should be taken as up- 
per limits, because of the heating of the donor star by 
the compact object. 

4. We calculated average absorption line profiles for the 
strong and weak lines separately, each line using 8 in- 
dividual lines. The position of the fine centers of the 
average absorption line agree within our CCF results. 
Prominent emission components are observed in the 
strong lines, indicating that the heating effects are im- 
portant for a proper interpretation. 

5. We construct a simple model where we take into ac- 
count the UV heating effects on the donor star surface 
and on its envelope, and where also consider the emis- 
sion from the gas stream. The model reproduces well 
the emission components and absorption Unes in the av- 
erage line profiles, both in intensity and radial velocity 
variations. These results indicate that the heating could 
have a significant impact on the estimate of the real ra- 
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dial velocity of the donor star, which may be as low as 
= 40 ± 5 km s"' . We then estimate the masses of the 
components as Mq = 10.4+f J Mq and Mx = 2.5!^;^ Mq . 

6. The final constraints for the compact object mass are 
1.9 Mq < Mx < 4.9 Mq, where the lower and up- 
per limits are inferred from the modeling of the aver- 
age absorption line profiles and from the CCF analy- 
sis, respectively. We conclude that the compact object 
in SS 433 is most likely a low mass black hole, al- 
though the possibility of a massive neutron star cannot 
be firmly excluded. 
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Fig. 1. — Optical photometry of SS 433 during tlie Subaru observations. Crosses show data points interpolated to the time of the spectral observations. 
Variations larger than 0.01 in V and 0.02 in B represent real photometric activity in SS 433. 
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Fig. 2. — Optical spectra of SS 433 observed with Subaru FOCAS covering the 3750-5250 A range. From top to bottom, the curves represent the spectra taken 
on the fourth ((/> = 0.262), third (0 = 0.1 15), first ((/> = 0.956), and second night (0 = 0.034). 
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Fig. 3. — "Normalized" spectrum of Region 1 taken during the first night of the Subaru observations. It contains several strong absorption Hues with emission 
components (cf. the highly rectified spectrum of Region 1 in Figure|5). The strong absorption line at A4501 is a diffuse interstellar band (DIB) blended with a 
Ti II line. 
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Fig. 4. — "Normalized" spectrum of Regions 2 and 3 taken on the first night of the Subaru observations. 
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Fig. 5. — Rest frame spectra of SS 433 (highly rectified) and HD 9233 (normalized) in Region 1 obtained with Subaru FOCAS. The SS 433 spectra are shifted 
by offsets of 0.4, 0.3, 0.2, and 0.1 for the first, second, third, and fourth night, respectively. Numbers on the right side indicate the corresponding orbital phases 
(p. The flux level of HD 9233 spectrum is reduced by a factor of 0.36. The absorption feature near 4501.79 A is of interstellar origin. 
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Fig. 6. — Same as Figure|3] but for Regions 2 and 3. The absoiption features near 4762.61 A, 4780.02 A, and 4963.88 A are of interstellar origin. 
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Fig. 7. — Highly rectified spectra of SS 433 obtained the Gemini GMOS in Region 2 and 3, together with the normalized spectrum of HD 9233 obtained with 
Subaru FOCAS. The absolution features near 4762.61 A, 4780.02 A, and 4963.88 A are of interstellar origin. 
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Fig. 8. — Radial velocity curve of the donor star in SS 433 obtained by the CCF analysis of the Subaru data, squares indicate the results from Region 1, filled 
circles from Region 2, and triangles from Region 3. 
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Fig. 9. — Radial velocity curve of the donor star in SS 433 obtained by the CCF analysis of the Gemini data. FiUed circles and triangles indicate the results 
from Region 2 and from Region 3, respectively. 
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Fig. 10. — Radial velocity curve of the donor star in SS 433 complied from both Subaru and Gemini results of the CCF analysis in Region 2. Filled circles 
and triangles represent Subaru and Gemini results, respectively. The best-fit Keplerian solution is displayed by the solid curve. Dashed lines correspond to the 
la error. 
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Fig. 11. — Observed average absorption-line profiles for the "strong" lines with emission components (right) and its best-fit model (left). From the first to the 
fourth Subaru night, the profiles are denoted by solid, dashed, dash-dotted and dotted Unes. 




Fig. 12. — Observed average absorption-line profiles for the "weak" Unes without notable emission components (right) and its best-fit model (left). From the 
first to the fourth Subaru night, the profiles are denoted by solid, dashed, dash-dotted and dotted lines. 
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Fig. 13. — Sketch of the binary system SS 433 (not to scale). The system consists of a mass donor star with an extended gas envelope, a UV heating source 
(surrounding the compact object), and a gas stream. The orbital motion and rotation of the donor are shown. We distinguish between three different regions of 
the donor surface, (1) the non-heated region from which we observe absorption lines only, (n) the heated region from which we observe emission lines only, and 
(111) the overheated region which does not emit any spectral lines. The extended envelope (i.e., the wind) produces emission lines only in Region 11 and in those 
parts of Region I that are exposed to the UV source. The orbital phases seen by the observer in the four Subaru nights are indicated. 
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Fig. 14. — Absorption-line radial velocity curves derived from the observed average line profiles (fiUed symbols) and those measured from the best-fit models 
using the same method as for the observed ones (open symbols). Circles represent the strong lines, and squares the weak lines, respectively. Kq = 40 km s~' and 
Kx = 160 km s~' are assumed. The heating model well reproduces the observed radial velocities. 
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Fig. 15. — Line profiles of He II A4686. Numbers shown in the upper right comer indicate the orbital phases. 
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Fig. 16. — Radial velocity curve of the compact object measured from the He II line. The horizontal line corresponds to the systemic velocity of ^0=59. 2 
km s~' as determined from the CCF analysis for the donor star. The dotted, sinusoidal curve represents the best-fit obtained using the data of the fourth night 
(Kx=159±7kms-i). 
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TABLE 1 

Log of Observations 



Date HID-2,450,000 


Total exposure 


Frames 




(Mid time) 


(sec) 




Subaru 








2007 Oct 06 


4379.76 


6000 


5 


2007 Oct 07 


4380.78 


8400 


7 


2007 Oct 08 


4381.84 


5484 


5 


2007 Oct 10 


4383.77 


8400 


7 


Gemini 








2006 Jun 07 


3893.99 


5400 


3 


2006 Jun 08 


3895.04 


9000 


5 


2006 Jun 09 


3895.96 


16200 


9 


2006 Jun 10 


3897.02 


12600 


7 


2006 Jun 11 


3898.01 


9000 


5 


2006 Jun 12 


3899.03 


7200 


4 


2006 Jun 13 


3900.03 


7200 


4 




TABLE 2 




Radial Velocity of the Donor Star Determined 


FROM THE CCF (Region 2) Together with 


PRECESSIONAL AND ORBITAL PHASES 




Date 


Observatory 


ijj (j> 


Vr 


(HJD-2,450,000) 






(kms"') 


4379.76 


Subaru 


0.018 0.956 


52±7 


4380.78 


Subaru 


0.024 0.034 


65 ±7 


4381.84 


Subaru 


0.031 0.115 


84±7 


4383.77 


Subaru 


0.043 0.262 


119±7 


3893.99 


Gemini 


0.022 0.838 


9±8 


3895.04 


Gemini 


0.029 0.918 


39±8 


3895.96 


Gemini 


0.034 0.988 


47±8 


3897.02 


Gemini 


0.041 0.069 


73±8 


3898.01 


Gemini 


0.047 0.145 


106±9 


3899.03 


Gemini 


0.053 0.223 


132±9 


3900.03 


Gemini 


0.060 0.300 


128±8 


Note. — The errors are la standard errors from the CCF 


analysis. 










TABLE 3 




Radial Velocity of the He II line 




Date 


Observatory 




Vr 


(HJD-2,450,000) 






(kms-l) 


4379.76 


Subaru 


0.018 0.956 


178±2 


4380.78 


Subaru 


0.024 0.034 


120±1 


4381.84 


Subaru 


0.031 0.115 


32±10 


4383.77 


Subaru 


0.043 0.262 


-96±7 



Note. — 



The errors are estimated systematic errors (see text). 



